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ABSTRACT 



o 



Aims. A nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants (SNRs) is employed to investigate the properties of 
SNR RX J 17 13.7-3946. 

Methods. Observations of the nonthermal radio and X-ray emission spectra as well as the H.E.S.S. measurements of the very high energy 
7-ray emission are used to constrain the astronomical and CR acceleration parameters of the system. It is argued that RX J 1 7 1 3. 7-3946 is a 
core collapse supernova (SN) of type Il/Ib with a massive progenitor, has an age of ~ 1600 yr and is at a distance of w 1 kpc. It is in addition 
assumed that the CR injection/acceleration takes place uniformly across the shock surface for this kind of core collapse SNR. 
Results. The theory gives a consistent description for all the existing observational data, including the nondetection of thermal X-rays and the 
spatial correlation of the X-ray and 7-ray emission in the remnant. Specifically it is shown that an efficient production of nuclear CRs, leading 
to strong shock modification and a large downstream magnetic field strength « 140 ^G can reproduce in detail the observed synchrotron 
emission from radio to X-ray frequencies together with the 7-ray spectral characteristics as observed by the H.E.S.S. telescopes. 
Conclusions. The calculations are consistent with RX J1713. 7-3946 being an efficient source of nuclear cosmic rays. 
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.£h 1. Introduction 

v : 

^ ■ RX J17 13-3946 is a shell-type supernova remnant (SNR) lo- 
' cated in the Galactic plane that was di s cover ed in X-rays 
with ROSAT dPfeffermann & Aschenbachl 1 1996b. Subsequent 



studie s of th is SNR with t he A SCA satellite by iKovama et al 
d 19971) 



and 



Slaneetal.1 dl999h . and later with XMM by 
Cassam-Chenai' et al.l ( 120041) . have not positively identified any 
thermal X-ray emission and therefore tentatively concluded 
that the observable X-ray emission is entirely non-thermal. 

The radio emission on the other hand is quite weak. Only 
part of the remnant shell could be detected in radio syn- 
chrotr on emission up to no w, with a poorly known spectral 
2004). For the spatially integrated radio 
2006) estimated about twice the value 



synchrotron flux increases the magnetic field estimate (see be- 
low). 

RX J1713-3946 was also detected in very high en- 
ergy (VHE:> 100 GeV) 7-rays with the CANGAROO 
dMuraishi et allbOOOtlEnomoto et all|2002b and H.E.S.S. tele- 
scopes dAharonian et al,, 2004, 2006, 2007). Especially the lat- 
ter, very detailed observations show a clear shell structure at 
TeV energies which correlates well with the ASCA contours. 



In a first theoretical paper on this source dBerezhko & Volk , 

2006) we have investigated the acceleration of CR elec- 



form ( Lazendic et al 



flux Aharonian etal 



found by lLazendic et al.l d2004l) . Most recently lAcero et al 
d2009l) derived a total flux value of 22Jy < S < 26Jy 
at 1.4 GHz, still abou t two times higher than estimated by 



Aharonian et al. (2006). As far as the estimate of the magnetic 
field strength is concerned, this higher estimate of the radio 
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trons and protons in detail, using nonlinear kinetic theory 
dBerezhko et all 1 19961: iBerezhko & Volkl l200d) . Observations 
of the nonthermal radio and X-ray emission spectra as well 
as the H.E.S.S. measurements of the very high energy 7-ray 
emission were used to constrain the astronomical and the par- 
ticle acceleration parameters of the system. Under the assump- 
tion that RX J1713. 7-3946 is the remnant of a core collapse 
supernova (SN) of type Il/Ib with a massive progenitor, the 
theory gives indeed a consistent description for the existing 
data on the nonthermal emission. Specifically the VHE data 
from H.E.S.S. were shown to be best explained as hadronic 7- 
ray emission, where the magnetic field amplification strongly 
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depresses the inverse Compton and Bremsstrahlung fluxes. 
Subsequently (IBerezhko & Volkl 120081) . we have analyzed the 
spatial correlation between the synchrotron and the VHE emis- 
sion that has been used befor e as an argument for a leptonic 
origin of the VHE emission dKatz & Waxmanl 120081: Plaga, 
2008). It was argued that correlated density and magnetic field 
strength variations lead rather naturally to a spatial correla- 
tion of the hadronic VHE emission with the synch rotron emis 



2. Results 



sion. S imilar arguments were brought forward by Tanakaetal 
(2008J). Also the recent broad-band X-ray synchrotron mea- 



surements dUchivama et all 120071; iTakahashi etaUl2008l) with 
the Suzaku instrument were compared to the theoretical syn- 
chrotron spectra and found to be quite consistent with a 
hadr onic model. A pur ely leptonic model on the other hand 
(e.g. lPorter et al. , 2006), where magnetic field amplification is 
not ex pected, was shown to be inconsistent with these observa- 
tions. (lAcero et al.l (120091) have also considered a leptonic ori- 
gin of the 7-ray emission as a viable possibility. This has to be 
compared with the results of the present paper.) 

These theoretical considerations assumed that only over 
part of the shock surface injection of suprathermal nuclear ions 
occurs effectively, because over other parts the magnetic field is 
essentially parallel to the shock surface. In fairly regular SNRs 
like those of Type la explosions into a presumably uniform 
circumstellar environment, this is an important effect, as for 
instance the case of SN 1006 shows. For a blast wave propa- 
gating into the highly turbulent shell of a stellar wind bubble 
from a massive progenitor star on the other hand, the mag- 
netic flux tubes in the shell are presumably so slender and 
irregular that energetic particles can cross effusively into re- 
gions where suprathermal injection is depressed, yet acceler- 
ation is possible. As a result shock modification by the pres- 
sure gradient of the accelerating particles becomes possible 
everywhere o n the shock sur f ace. T his reasoning follows the 
discussion in Berezhko et al. I d2009l) with the conclusion that 
SNRs propagating into wind bubbles should have their forward 
shock modified everywhere. Technically this implies that the 
SNR bec omes spherically symmetric and the correction factor 
/ rc < 1 dVolk et all 12.0031) on the spherically symmetric solu- 
tion reaches unity. The ubiquitous shock modification has also 
implications for the thermal emission, because the thermal gas 
is then heated only in the subshock and the average gas density 
required for a given hadronic 7-ray emission is reduced. 



Contrary to the previous studies dBerezhko & Volki r2006. 



20081) of RX J1713.7-3946 we therefore adopt here the approx- 
imation of spherically symmetric CR in jection/acceleration 
with / rc w 1 (see also Zirakashvili, 2009). It is demonstrated 
that also in this approximation a consistent solution can also be 
found. In addition a rough estimate of the thermal X-ray emis- 
sion from this object shall be given. Even though the error in 
this estimate might be quite large, the nominal result is con- 
sistent with the pr esent non-detection of thermal X-rays (e.g. 
Acero et all 120091) . Finally, the observed correlation of the X- 
ray and the VHE 7-ray emission is discussed in a generalized 
form. 



As in the previous studies dBerezhko & Volk , 2006, 2008) a 
source distance of d = 1 kpc is adopted here. For the present 
angular source size of about 60 arcmin this implies a SNR ra- 
dius R s ss 10 pc. 

For the SNR age th e value t sn = 1612 yr is used, consistent 
with the hypothesis of IWang et all dl997b that RX J1713-3946 
is the remnant of the AD 393 "guest star". 

It had been argued (Berezhk o & Volkl 120061) that the ex- 
isting data are only consistent with RX J1713. 7-3946 being 
a core collapse SN of type Il/Ib. Those progenitor stars of 
core collapse SNe that emit intense winds are massive main- 
sequence stars with initial masses Afj ^ 15Mg. During their 
evolution in a surrounding uniform ISM of gas number density 



= m„N-\ 



the wind termination shock creates a hot bub- 



ble of shocked wind material, surrounded by a turbulent shell 
of shocked ISM behind an outer forward shock that communi- 
cates the internal overpressure to the environment. 

The gas number density distribution N g (r) — g(r)/m p is 
assumed to have the form 



N g = 0.008 + 0.24[r/(10pc)] 



12 cm" 3 



(1) 



This density profile corresponds to a gas number density N g = 
0.25 cm -3 at the current SN shock positio n, which is a factor 
~ 5 lower than in lBerezhko & Volkl d2006l) . 

Under these conditions a SN explosion energy E sn = 
1.3 x 10 51 erg and an ejecta mass M e j = 3M Q lead to a good 
fit for the observed SNR properties. To determine the explosion 
energy directly from the data the value of the shock speed V s 
would have been needed. However, this quantity is not mea- 
sured. The ejecta energy E C j is still about 0.25-E sn . Together 
with the kinetic energy of the shocked gas E^ w 0.22? sn this 
makes up about half of the total energy E sn (see Fig.l.). 

The spectral fit, also for the 7-ray data, can be achieved 
with the values r\ = 5 x 10~ 4 of the proton injection rate and 
Bd = 142 /iG for the downstream magnetic field strength. 
Such values of £?d are significantly higher than typical ambi- 
ent magnetic fields, also for the considered bubble wall in the 
dense environment (the gas density at the shock is still a factor 
~ 100 lower than in that environment which implies a strongly 
decompressed field upstream of the SNR shock compared to 
the value in the ambient medium). The high -value must be 
attributed to field amplification at the shock front due to the 
strong wave produ ction by the acceleration of C Rs far into the 
nonlinear regime dLucek & Bell 120001 belli |2004 Q 

It is assumed that also electrons are injected into the ac- 
celeration process at the shock front. Since the details of the 
electron injection process are poorly known, the electron injec- 
tion rate is chosen such that the electron to proton ratio K ep 
(which is defined as the ratio of their distribution functions at 
all rigidities where the protons are already relativistic and the 
electrons have not been yet cooled radiatively) is a constant 
to be determined from the synchrotron observations. Clearly, 



1 For direct observational eviden ce regarding magn e tic field ampli - 
fication in RX J1713.7-3 946 see IVolk etal.1 d2005h ; iBalled d2006h ; 
lBerezhko& Volkl j2006h . 
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Fig. 1. (a) Shock radius R s , shock speed V s , and preshock gas 
density 7Y g ; (b) total shock (er) and subshock (cr s ) compression 
ratios; (c) ejecta (E e y) gas kinetic (Eg0 and CR (E c ) energies 
as a function of time, normalised to the total mechanical SNR 
energy E sn . The vertical dotted line marks the current evolu- 
tionary epoch. 



from the point of view of injection/acceleration theory, K cp , 
together with Bo = B^/cr and r\, must be treated as theoreti- 
cally not very well constrained parameters to be quantitatively 
determined by comparison with the available synchrotron and 
7-ray data. Here a is the overall shock compression ratio. 

In the present case the following parameter values are ob- 
tained by iteration: 77 = 5 x 10~ 4 , Bq = 25 /iG, K cp = 
1.3 x 10~ 3 . The gas dynamic variables at the present epoch 
are a = 5.7, er s = 3.3, and V s w 2200 km sec -1 , where a s 
and V s denote the subshock compression ratio and the overall 
shock velocity, respectively (Fig. 1). 



The corresponding solutions of the dynamic equations at 
each instant of time yield the CR spectrum and the spatial dis- 
tributions of CRs and thermal gas, and therefore also the ex- 
pected fluxes of nonthermal emission produced by the acceler- 
ated CRs. 

The overall broadband spectral energy distribution (SED), 
expected at the current evolutionary epoch, is displayed in 
Fig. 2, together with the experimental data from A TCA at radio 



and Suzaku (Uchivama et al. 
data dAharonian et alll2007h . 



wavele ngths, as estimated for the full remnant b y Acero et al 
d2009h . the X-ray data from ASCA dAharonian et all 12006)) 



2007), and the H.E.S.S. 7-ray 



The comparison wit h our earlier theoretical spectrum 



dBerezhko & Volk , 2008 ) shows that also a model with full 
spherical symmetry can successfully describe the spatially in- 
tegrated nonthermal emission properties of this object. 

At energies e 7 < 1 TeV the theoretical 7-ray spectrum is 
as hard as dF 7 /de 7 oc e -18 , whereas for e 7 > 10 TeV it has a 
smooth cutoff. It should be noted that the 7-ray cutoff energy 
£ max 0.1cp max is sensitive to the magnetic field strength B^, 
since the proton cu t off m omentum has a dependence p max oc 
R s V s B d dBerezhko[[l996l) . 

The theoretical results presented here are fully consistent 
with a dominantly hadronic origin of the observed TeV 7-rays 
. At the present epoch the SNR has already converted ss 35% 
of E* n into accelerated nuclei. 



2. 1. Estimate of the thermal X-ray emission 



Like for the case of RX J0852.0-4622 (Berezhkoetafll2009) 



in addition the expected flux of thermal X-rays is estimated 
here. This is only possible in a very approximate way, even if 
we shall disregard the ejecta emission with the argument that 
the ejected mass in the present phase is small compared to the 
swept-up mass. The reason is the explosion into the wind bub- 
ble, creating a rather different thermodynamical structure from 
that of a classical Sedov solution for a SN explosion into a 
uniform medium, exclusively until now considered in the lit- 
erature. Yet the remnant of RX J1713. 7-3946 is well past the 
sweep-up phase and has entered a quasi-Sedov phase in the 
stellar wind shell, i.e. a roughly self-similar evolutionary phase, 
modified by strong particle acceleration relative to a purely gas 
dynamic evolution. Other estimates of the thermal emission 
from R X J 17 13.7-3946 have been made by iKatz & Waxmanl 
( 2008) and bruryet all d"2008 ). effectively assuming an explo- 
sion into a uniform medium. 

The approximation adopted here is the following. The 
present bubble case is compared with a SNR in a uniform 
medium in the classical Sedov phase without any CR acceler- 
ation, making four assumptions (i) the total hydrodynamic ex- 
plosion energy is the same in both cases (ii) the shock velocity 
is the same (iii) the present gas density upstream of the shock is 
the same, and (iv) the t wo objects are at the sa me distance of 1 
kpc. Then the results of Hamilton et al.l d 19831) for the thermal 
X-ray flux from a classical Sedov SNR are used, employing the 
emission measure of the bubble remnant instead of that of the 
classical Sedov remnant with the same four parameters above. 
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log e , eV 



Fig. 2. Spatially integrated, overall nonthermal spectral energy distribution of RX J 17 13.7-3946 The solid curve at energies above 
10 7 eV corresponds to 7r°-decay 7-ray emission, whereas the dashed curve indicates the Inverse Compton (IC) emission. The 
dotted line corresponds to the test particle limit which implies insignificant proton acceleration and m agnetic field amplification 



(seelBerezhko & VolklEooI for the details). T he ATCA radio data, as derived by lAcero et ail J2009t) . th e ASCA X-ray data (cf . 
1 120061) . 



Aharonian et al. 



the Suzaku X-ray data dUchivama et all 120071) . and the 2006 HESS 7-ray data dAharonian et all 2007) 



are also shown. The EGRET upper limit for the RX J1713. 7-3946 position dAharonian et all 2006) is included as well. 



This means that the X-ray emissivity of the remnant is reduced 
by the ratio i? cm = EM^ / EM§ of the emission measure EM\, 
for our bubble solution to the emission measure for the classical 
Sedov solution EM s, which corresponds to a uniform ambient 



gas density N g (R s ) (Berezhko et al 
The relation 

T sub _ [(a s -l)( 7 +l) + 2]( 7 + l) 2 



2009). 



F 



47(7 - l)er 2 



(2) 



is used, with 7 = 5/3, to determine the subshock temperature 
T su b, relevant for the thermal emission flux for the modified 
shock. T s = 10 7 (V;/(839km/s)) 2 K w 6.9 x 10 7 K denotes the 
gas temperature of an unmodified shock with V s « 2200 km/s. 

Using the differential thermal X-ray model spectra 
dF/de w 3 x 10" 5 photons/(keV cm 2 s) from Hamilton et al 



( 1983J) (for T e ^ T-, see their Fig.2) for their 17 = N&E m = 



10 erg cm" b with account of the scaling dF/ de oc 77 with our 
model numbers 77 w 8 x 10 49 erg cm" 6 , multiplying it by the 
factor 6 2 [E sn / (lO^erg)]" 1 / 2 , as required, where 8 = 54.5 ar- 
cmin for the angular size of the classical Sedov remnant corre- 
sponding to RX J1713. 7-3946 and multiplying it also by the 
factor i? om = 0.46, results in a thermal spectral energy density 
e 2 dF/de « 288 eV cm" 2 sec" 1 for e = 1 keV. 

This result has not yet taken into account that the actual 
postshock gas temperature T su b is by a factor rs 0.41 lower 
than the temperature Ts used i n the above estimate. As dis- 
cussed in Ber ezhko et all (120091) this roughly corresponds to a 
reduction of the parameter value 77 rs 3.3 x 10 49 erg cm" 6 and 
thus to a reduction of the thermal flux at 1 keV to e 2 dF/de ~ 
118 eV cm" 2 sec -1 . 



These numbers must be compared with the observed non- 
thermal X-ray energy flux at 1 keV (see Fig.l) e 2 dF/de rs 
200 eV cm" 2 sec" 1 . Therefore, at 1 keV, the thermal energy 
flux comes out to be smaller than the nonthermal flux by a fac- 
tor of about 2. This does not contradict the non-detection of 
thermal X-rays from RX J 17 13.7-3946 . 

2.2. Spatial correlation of the nonthermal X-ray 
fluxes with the j-ray fluxes 



Note that Acero et al. (2009), who find a positive spatial non- 



linear correlation Fx oc F" with a = 2.41 between the flux F 1 
of 7-rays with energies e 7 = 1 — 10 TeV and the flux of non- 
thermal X-rays Fx with energies e 1 = 1 — 10 keV, concluded 
that such a correlation is in favor of a leptonic nature of the 
TeV-emission. Such a positive correlation is naturally expected 
within the present approach. 



2.2.1 . Correlation of the radial profiles 

In terms of the spherically symmetric model, discussed here, 
the correlation between X-ray and 7-ray fluxes means that the 
projected radial profiles Jx (p) and J 1 (p) of these emissions as 
functions of projected radial distance p have a similar shape. At 
first sight this is a contradiction to this model: as can be seen 
from Figs. 3 and 4, the calculated profiles are significantly dif- 
ferent. Both profiles have a sharp peak just behind the shock 
front, but the X-ray profile is by a factor of about 3 thinner 
than the 7-ray profile. However, when smoothed with the point 
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p/R s 

Fig. 3. The X-ray emissivities for the energy e 7 = 1 — 10 keV 
as a function of projected, normalized radial distance p/i?, s .The 
calculated radial profile is represented by the solid line. The 
dashed line represents the calculated profile convolved with the 
point spread function of Gaussian width 0.05°. 




p/R s 

Fig. 4. The 7-ray emissivities for the 7-ray energy e 7 = 
1 — 10 TeV as a function of projected, normalized radial dis- 
tance p/R s . The calculated radial profile is represented by the 
solid line. The dashed line represents the calculated profile con- 
volved with the point spread function of Gaussian width 0.05°. 



spread function of width do = 0.05°, corresponding to the 
H.E. S.S. angular resolu tion, these profiles become very simi- 



lar. dAcero et all 120091) compared XMM measurements of the 



X-ray emission smoothed with a point spread function of width 
a = 0.05°. It can therefore be concluded that the similarity 
of the X-ray and the 7-ray radial profiles is a trivial conse- 
quence of the smoothing of two differently sharp peaks by a 
point spread function whose width is considerably larger than 
the widths of either one of the emission peaks. 

The experimental ly observed correlation according to Fig. 9 



of Acero et al. I J2009I) is rather close to the relation Fx = F 7 



for all the brighter SNR regions. Some deviation from this re- 
lation is observed for six of the outer and dimmer regions, for 
which Fx < F T This peculiarity can be also explained within 
the present model. The reason is that protons with an energy 
of about 100 TeV occupy a region, which is noticeably larger 
than the shock size R s , due to their high diffusive mobility. 
The 10 TeV 7-ray emission which these protons produce has 
an even larger radial extension since the gas density strongly 
increases with radial distance. Such an effect is absent for the 
electrons of the same energy which produce the nonthermal 
keV-emission: as a result of their strong energy losses they 
occupy a region thinner than 10 _2 i? s around the shock front. 
Therefore, as can be clearly seen from Figs. 3 and 4, the ratio 
Fy /Fx is expected to incre ase as a func t ion of radial distance 
for p/R s > 1, as the data of Acero et al. ( 20091) are interpreted 
here. 



2.2.2. Correlation of azimuthal variations 

The measured X-ray and 7-ray emissions also undergo rather 
closely correlated azimuthal variations across the remnant. 
Such an effect can not be described within the present model 
model, because it is spherically symmetric. One can neverthe- 
less attempt to interpret the correlated behaviour of the fluxes 
Fx and F 7 in different azimuthal sectors of the remnant. It is 
natural to assume that this variation of the remnant properties 
is due to variations of the ambient gas density g. According to 
the model proposed here the SNR's present evolutionary phase 
is intermediate between the sweep-up and a quasi-Sedov phase. 
Since the ejecta kinetic energy together with the kinetic energy 
of the swept-up gas contain about half of the explosion energy, 
the SN shock can be treated approximately as a piston-driven 
shock. The speed of such a shock depends only weakly upon 
the upstream gas density. Therefore can be approximated by 
a constant across the whole remnant. 

The local hadronic 7-ray emission then varies like F 7 oc 
qE c oc g, where the local CR energy content E c in each angu- 
lar sector is assumed to be the same everywhere, because the 
overall CR energy content E c k, 0.4Fsn has already reached 
the saturation level and therefore can not undergo significant 
variations. In addition, the highest energy CRs in the SNR have 
the highest diffusive mobility so that they tend to be distributed 
roughly uniformly across the remnant. 

As it is clear from Fig. 2, the synchrotron spectrum at high 
energies e 7 > 10 eV is already considerably influenced deter- 
mined by synchrotron losses of the emitting electrons. Since 
the energy of these electrons is rather rapidly and completely 
transformed into synchrotron emission, the flux Fx oc K cp gV^ 
is determined by the available energy flux pV^/2 into the 
shock, an essential part of which is transformed at the shock 
front into the energy flux of high energy CRs. High energy 
electrons therefore accumulate a fraction K cp of this flux. 
This flux hardly depends on the value of the interior mag- 
netic field in all those SNR regions where the field is suffi- 
ciently high, Fd > 140 pG. In this case the relation Fx oc F 7 
is expecte d for the br i ghtest part of remnant. This is indeed 
observed dAcero et all [2009). Since the energy density of the 
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amplified magnetic field is phenomeno logically known t o 
be proportional to gVf, with 2 < (3 < 3 dVolk et all 120051) . 
it is expected to be proportional to the local upstream gas 
density g. The magnetic field strength is therefore lower than 
the average value w 140 fiG within those parts of the rem- 
nant, where the gas density is lower than the average value. 
Synchrotron losses become much smaller in these regions and 
therefore the expected synchrotron X-ray flux is approximately 
Fx oc gK cp B^ 2 oc g 7 / 4 . For the dim part of the remnant, in 
which X-ray and 7-ray emission is considerably lower than on 

7 /4 

average, one then obtains Fx oc Fj . This is approximately 
consistent with the observational result. 

3. Summary 

The assumption that CR injection/acceleration takes place uni- 
formly across the whole SN shock surface is consistent with 
the existing data. 

The swept-up mass is so small in this case that, in a rough 
approximation, the estimated flux of thermal X-rays at 1 keV is 
lower than the nonthermal flux, consistent with the nondetec- 
tion of thermal X-ray emission until now. 

It is concluded that the present observational knowledge of 
SNR RX J1713. 7-3946 can be interpreted by a source which 
ultimately converts more than 35% of the mechanical explosion 
energy into nuclear CRs. Also the observed high energy 7-ray 
emission of SNR RX J1713. 7-3946 turns out to be primarily 
of hadronic origin. 
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